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We present a novel analysis of diffractive DIS data, in which the input parton distributions of the Pomeron are 
parameterised using perturbative QCD expressions. In addition to the usual two-gluon model for the perturbative 
Pomeron, we allow for the possibility that it may be made from two sea quarks. In particular, we treat individually 
the components of the Pomeron of different size. This property allows the absorptive corrections to the inclusive 
DIS structure function F2 to be calculated using the AGK cutting rules. The absorptive effects are found to 
enhance the size of the gluon distribution of the proton at small x. 



1. DIFFRACTIVE DIS 

A notable feature of deep-inelastic scattering is 
the existence of diffractive events, 7*p — > Xp, in 
which the slightly deflected proton and the cluster 
X of outgoing hadrons are well-separated in ra- 
pidity. The large rapidity gap is believed to be as- 
sociated with Pomeron exchange. The diffractive 
events make up an appreciable fraction (« 10%) 
of all (inclusive) deep-inelastic events, j*p — * X. 
We will refer to the diffractive and inclusive pro- 
cesses as DDIS and DIS respectively. 

The diffractive structure function, 
F2^^\xip, (3,Q'^), has recently been measured 
to high precision by the ZEUS [H and HI (2] 
Collaborations at HERA 3|. 

Both HERA experiments have analysed their 
DDIS data using the 'resolved Pomeron' model 
where the xip dependence of F^''^'' is assumed 
to factorise into a 'Pomeron flux factor' and the 
Pomeron is treated as having parton distribution 
functions (PDFs) just like a hadron. Since the 
Pomeron is not a particle, this model must be 
taken on a purely phenomenological basis. 

An alternative description of DDIS is provided 
by interpreting Pomeron exchange as two-gluon 
exchange, such as in the BEKW model 5 . How- 
ever, the BEKW parameterisation is rather far 
from the original perturbative QCD (pQCD) cal- 
culations of Wiisthoff |6j . 



2. NEW pQCD APPROACH TO DDIS 

In pQCD, it is known that the Pomeron singu- 
larity is not an isolated pole, but a branch cut, in 
the complex angular momentum plane [Jj. The 
pole singularity corresponds to a single particle, 
whereas a branch cut may be regarded as a con- 
tinuum series of poles. That is, the Pomeron wave 
function consists of a continuous number of com- 
ponents. Each component i has its own size, l//ii. 
The DGLAP evolution of a component should 
start from its own scale /ii, provided that fii is 
large enough for the perturbative evolution to be 
valid. Therefore, the expression for the diffrac- 
tive structure function F^^^"^ contains an integral 
over the Pomeron size, or rather over the scale fi: 




xfJp{xIP■,^i')F^{(3,Q^■,^l^). (1) 

Here, the subscript P on F^p^^ is to indicate 
that this is the perturbative contribution with 
fi > Qo ^ 1 GeV. The notation IP = G, S, GS 
denotes that the perturbative Pomeron is repre- 
sented by two i-channel gluons, two i-channel sea 
quarks, or the interference between these, respec- 
tively; see Fig.n The necessity for the 'two-quark 
Pomeron' in addition to the two-gluon Pomeron 
is due to the valence-like shape of the gluon distri- 
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Figure 1. New pQCD approach to DDIS. Each 
component of the perturbative Pomeron of size 
is represented by either (a) two i-channel 
gluons in a colour singlet or (b) sea quark- 
antiquark exchange. The perturbative Pomeron 
flux factors fip{xjp; fj?) are given in terms of the 
gluon and sea quark distributions of the proton, 
g{xip,fi^) and S{xjp,^'^). The Pomeron struc- 
ture function (/3, Q^; /i^) is evaluated at NLO 
from the quark singlet, TF{j3, Q^; ^^), and gluon, 
g^iP, Q^; fi^), distributions of the Pomeron. 



bution of the proton at low scales. The (3 depen- 
dence of the input Pomeron PDFs, S^(/3, y-'^^H^) 
and g^{(3,fi^;fi^), is obtained from the lowest- 
order Feynman diagrams, with the normalisations 
allowed to go free in fits to the DDIS data to ac- 
count for higher-order QCD corrections. 

In addition to the perturbative contribution flj 

to F^^'^\ we also include a non-perturbative (NP) 
contribution, a twist-four contribution, and a sec- 
ondary Reggeon (IR) contribution. For more de- 
tails, see 

3. ABSORPTIVE CORRECTIONS TO Fa 

The total proton structure function, 
F2'^^^{xb,Q^), as measured by experiment, can 
be approximately written as a sum of the one- 
Pomeron (DGLAP) contribution and absorptive 
corrections due to two-Pomeron exchange; see 
Fig. 13 Apphcation of the AGK cutting rules 
to the IP (X> IP contribution gives 
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Figure 2. Absorptive corrections to F2 due to the 
IP® IP contribution, and application of the AGK 
cutting rules [5]. 



where F^ {xb,Q'^; ^'^) is the contribution to the 
diffractive structure function F^^^^ (integrated 
over xjp) which originates from a perturbative 
component of the Pomeron of size l//i. We can 
thus separate the screening corrections coming 
from low fj, < Qo, which are included in the in- 
put parameterisations taken at a scale Qq, from 
the absorptive effects at small distances (/i > Qo) 
which occur during the DGLAP evolution in the 
analysis of DIS data. For more details, see 

4. SIMULTANEOUS F2 + ^2^'^' FIT 

The 'simultaneous' fit of DIS and DDIS data 
proceeds as follows: 

(i) Start by fitting ZEUS [TTl and HI |T2 F2 
data (279 points) with no absorptive correc- 
tions, similar to the MRST2001 NLO anal- 
ysis [Tl^ . 

(ii) Fit ZEUS mi and HI [HI F^^^^ data (404 
points) using g{xjp,fi'^) and S{xip,ix^) from 
the previous F2 fit. 

(ui) Fit F2°°L^P = F2^ata-t- 1 AFa^'''^ | , with AFf''" 
from the previous F^^^"^ fit. 

(iv) Go to (ii). 
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x„ = 0.00005 to 0.00032 



xg(x,Q"=l GeV") = x" ■ (1-x) ■ (1 + E^x ■ ) - A x" (1-x) 
xS(x,Q^=l GeV') = A a''' (1-x)" (1 + E x°') 
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Figure 3. F2 data at small Xb before and 

after absorptive corrections have been included. 
Only data points included in the fits, shown by 
the solid and dashed lines, are plotted. 



In practice, only a few iterations of steps (ii) and 
(iii) are needed for convergence. 

In Fig. Owe show the F2 data at the smallest 
Xb values, before and after the absorptive correc- 
tions have been applied. The predictions of the 
corresponding fits, shown by the solid and dashed 
curves, are also plotted. In Fig. ^ we show the 
input PDFs at Qo = 1 GeV^, and also the gluon 
distribution when evolved up to = 4 GeV^. 
The dashed curves in Fig. 0] show that the in- 
clusion of absorptive effects yield an input gluon 
distribution which is much less negative, whereas 
the input sea quark distribution is largely unaf- 
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Figure 4. The gluon and sea quark distributions 
obtained from a NLO DGLAP fit to F2, before 
and after absorptive corrections have been in- 
cluded. 



fected. In fact, repeating the fits without the neg- 
ative term in the input gluon distribution gives 
a description of the F2 data which is almost as 
good (x^ /d.o.f. = 1.11, compared to 1.09 for a 
negative input gluon). By contrast, without any 
absorptive corrections, the fit to F2 is much worse 
without the negative term (x^ /d.o.f. = 1.57, com- 
pared to 1.15 for a negative input gluon). We 
conclude that absorptive corrections remove the 
need for a negative gluon distribution at Qq — 1 
GeV2. 

The corresponding description of the DDIS 
data, with a x^ /d.o.f. = 1.14, is shown in Fig. [S] 

5. DIFFRACTIVE PDFS 

We can extract diffractive PDFs (DPDFs) from 
the leading-twist contributions to the F^''^'' fit 
shown in Fig.jS] The diffractive quark singlet and 
gluon distributions are plotted in Fig. El and com- 
pared to those obtained by HI I2il5j . Note that 
the diffractive quark singlet distribution obtained 
by HI has a slightly steeper dependence, and 
hence HI obtain a larger diffractive gluon distri- 
bution. In addition, the smaller value of as{M'^) 
used by HI also enlarges their gluon density. It 




Figure 5. Description of combined preliminary ZEUS ^31 and HI jJBl DDIS data. 
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Figure 6. DifFractive PDFs, including a break- 
down of the different components, compared to 
those obtained by HI )2I15| . 



7. CONCLUSIONS 

In summary, we have achieved a good simul- 
taneous description of all the DDIS and small- 
Xb inclusive DIS data, in which the absorptive 
corrections in the description of the latter data 
have been identified and incorporated. In this 
way a more physical input gluon distribution 
at = 1 GeV^ has been obtained, which no 
longer needs to be negative at small x. However, 
there remains an outstanding dilemma in small-x 
DIS. Either, contrary to expectations, the non- 
perturbative Pomcron does not couple to gluons, 
or DGLAP evolution is frozen at low Q^, per- 
haps by power corrections. Note, however, that 
in both scenarios we still have the puzzle that the 
secondary Reggeon couples more to gluons than 
to sea quarks. 



has been demonstrated by HI that their DPDFs 
[211 5 j can be used to describe final state observ- 
ables in DDIS, namely dijet and D* meson pro- 
duction cross sections ^^l- Before our DPDFs 
can be taken seriously we need to demonstrate 
this also. 

6. POMERON-LIKE SEA QUARKS BUT 
VALENCE-LIKE GLUONS? 

The inclusion of absorptive corrections has not 
removed a long-standing puzzle concerning the 
behaviour of PDFs at small x and low scales. 
We still have a valence-like gluon distribution, 
whereas the sea quark distribution increases with 
decreasing x. That is, since the HERA F2 data 
have become available, we have had a 'Pomeron- 
like' sea quark distribution. According to Reggc 
theory, the high energy (small x) behaviour of 
both gluons and sea quarks is controlled by the 
same rightmost singularity in the complex angu- 
lar momentum plane, and so we would expect 
\g = Xs, where the input xg ^ x'^o , xS ^ x~^^ . 
We have studied several possibilities of obtaining 
a satisfactory fit with this equality imposed. The 
only modification which appears consistent with 
the data is the inclusion of power-like corrections, 
which may slow down the DGLAP evolution at 
low below about 1 or 2 GeV^. 
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